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SUMMARY 



A method of calculating the stresses and torsional 
stiffnesses of thin shells with interior webs is summa- 
rized. Comparisons between experimental and calculated 
results are given for 3 duralumin beams, 5 stainless steel 
beams, and 2 duralumin wings. It is concluded that if the 
theoretical stiffness is multiplied by a correction factor 
of 0.9, experimental values may be expected to check cal- 
culated values within about 10 percent. 



INTRODUCTION 



It is well known that the application of ordinary en- 
gineering formulas to thin sheet-metal structures may lead 
to serious errors in some cages. Eor certain types of cal- 
culations, however, standard formulas will give a reasona- 
ble degree of accuracy if experimentally determined correc- 
tion factors are applied. In the present paper U.S. Army 
Air Corps data on the tests of wings of the central box 
type are analyzed and the correction factor for initial 
torsional stiffness is derived. Since the methods of cal- 
culating torsion tubes with interior walls are not very 
widely known, the' first part of the paper gives a summary 
of a convenient method elucidated by a numerical example. 
The second part of the paper discusses the most important 
features of. the test objects and results. 
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LIST OF SYH30LS 

T, torque, in.-lb. 

f g , shearing stress, I'd. per sq.in. 

t , thicknes s , in. 

ds, differential element of perimeter, in. 

A, area bounded by median line, sq.in. 

J, torsion constant, analogous to moment of iner- 
tia in bending, in.* 

G, modulus of shear, lb. per sq..in» 

G J , torsional stiffness, lb. -in. 8 

9, angle of twist, radians (per unit length). ■ 

T| > factor of effectiveness in torsional stiffness, 

I. THE0R3T I OAL FORMULAS 



General remarks .- The fundamental torsion formulas for 
thin shells such as shown in figure 1 are 

f = (1) 
s - 2At V ' 

4A S 

J = TIE <*> 

the integral being taken around the perimeter of the pro- 
file. The derivation and the assumptions underlying these 
formulas nay be found in any good textbook on strength of 
materials, such as reference 1. 

The corresponding formulas for shells with longitudi- 
nal interior webs (fig. 2) are not so well known as the 
fundamental formulas. Their derivation by different meth- 
ods may be found in references 2 to 5, of which 2 and 5 
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are nost readily available. The present paper will con- 
fine itself to showing how the torsional stresses and 
stiffnesses of such shells may be calculated. She method 
presented here is that given in reference 4 and was choson 
"because it permits the shortest possible exposition of the 
subject. The method givon in references 2 and 5 is very 
similar; reference 3 giveSj a method based on the membrane-* 
analogy (reference l) . The most elaborate discussion of 
the physical foundations of the theory is given in refer-"" 
ence 5. 

M ethod o f calculation .- The interior webs will divide 
the cross section, of the shell into a number of separate 
cells. lumbers from one to n are assigned to these 
cells; the number zero is assigned to the space- outside of 
the shell. The wall and anything pertaining to it between 
the two cells i and j are designated by the subscript 
ij*. 5"o r each cell i is now computed the area Aj_ and 
for each wall ij the line integral ' " ~"" * 



a 



An auxiliary function T is now introduced for each cell 
i except cell zero, and the -equation 

i=i _ .... 

and n equations of the form . 

~ S (»j .- ^i) aij + 2&9 =0 (5) 

may be written down. The summation in equation (5) ex- 
tends over - all walls bounding the cell i, i.e., for j 
are substituted in turn all numbors between "iTero and n 
that designate the cells surrounding cell i. 

The system of equations (4) and (5) is solved for the 
unknowns G-8 and The. sheading stress in the wall 

ij is then given by 
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The direction ij is positive if a point traveling along 
the wall keeps cell i on its left. In figure 3 the pos- 
itive directions are indicated "by arrows. It is evident 
that 



f = -f „ 

s i«3 s Ji 



The torsion constant follows from 

J = 



T 

G-9 



(?) 



(8) 



The equations (4) and (5) applied to the three-cell 
system of figure 2 are 



T = 2 + A 3 1*3 -f-A 3 T 3 ] 



(4a) 



J--[-a 30 F 3 + a 31 (2' 1 - 3F 3 ) + a s3 (F 3 " %)] + 2&9 = 0 j-< 5a > 

j£-[-*so*a + a 3S t3P a - I 3 )] = 0 

For a shell with only one interior web (fig. 3), the 
values of J\ and J can be obtained directly from 



IS 



1 

2 



1 

2 



4 



l 2 0 



A 1 + a 12 A 



a s o l ®a a A a o l 



T 



a c l Aa 9-1 a A 



13 J 



a S0 A l 3 + & \3^ + a Ol A J 



a oi a i2 + a -i=i a =o + a^na 



13 a 30 



*-3 0 a O 1 



where A 
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For more than one interior web the formulas become too 
cumbersome, and it is more convenient to write down th e ~ 
equations (4) and (5) after the coefficients have been ob- 
tained and to solve these numerical equations. 

The method discussed can be used for analyzing a 
shell such as shown in figure 4 where the cover sheet is 
provided with closed longitudinal stiffeners. Co nsi der 
firsV the case in which tho top and bottom sheets areTbf 
the same thickness and all stiff oners are of the same size 
(fig. 4(a)). As indicated on "the figure by dotting all 
but one of the n stiffeners, only one stiffener need be 
considered because all functions will be the same for all 
stiffeners. Instead of obtaining a system of (n + l) 
equations of the type (5) , only the following three equa- 
t ions appear : ... _ . ... 

I = 2 [3?^ + n^Aj 



^-[-a 30 ? s + a 31 (F 1 - 1,)]+ 309 - ' 0 



> 



(10) 



where A 2 = A - nA 3 , A denoting the area bounded by the 
outer shell and A 3 the area under one stiffener; the 
line integrals a 30 and a ia are taken for one stiffener. 



In the more usual case where the top and bottom covers 
have different thicknesses and stiffeners, the equations 
become 



* = 2 [P X J^ + nTaAg + nEjAg] 



• I x ) + ma 13 (F 3 
% ) j + 2S9 = 0 



F, ) + 2G6=0 



(11) 



~-[-a 30 F 3 + a 31 (F 1 - 3? 3 )] + 2G9 = 0 



G 
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In these equations, the cell number- 1 refers to the main 
part of the shell, number 2 to the area under a top stif- 
fener (as in fig. 4(a)), and 3 to the area under a bottom 
stiffener; there are n top stiffeners and m bottom 
stif f ener s. 

When computing the line integral aj_ for a stiffener, 
it is conservative to use the full developed length and un- 
conservative to use the developed length botween rivets. 
If the average of the two lengths is used, the error in J 
should be very small bocause stiffeners contribute normal- 
ly only little to the total torsional stiffness. 

The equations for the case of figure 4(a) may be ob- 
tained in a slightly different manner that is, perhaps, a 
little more physically obvious. Since the stress condi- 
tions are equal and uniform along the top and bottom sur- 
faces, the stiffness of the shel* is obviously not changed 
if all the stiffeners are transferred to one side and com- 
bined into a continuous corrugated sheet (fig. 4(b)). All 
rivet rows except the outer two may now bo removed and the 
shell is reduced to a shell with a single Interior web for 
the purpose of computing the initial stiffness. Naturally, 
such a relocation of stiffeners would change the stiffness 
under large torques as well as the strength. 

Ap-pro xi mate method fo r stiffened, cover- .- The calcula- 
tion of the torsional stiffness by the method- described 
requires, in some cases, more time than is warranted con- 
sidering the probable magnitud-o of the discrepancies bo- 
tween calculations and tests and considering the purpose 
of tho calculation. In all cases where only an estimate 
of the torsional stiffness is required, and possibly in 
other cases, it will be sufficiently accurate to use the 
following approximate method for shells with stiffened 
cover . 

Find the centroid of a single stiffener (without skin 
attached to it) and join the centroids of all stiffeners 
on one cover sheet by a line. 

Along the line thus established, distribute the "ef- 
fective" matorial of__ all stiffeners. The effective mate- 
rial of one stiffener is the actual material multiplied by 

the ratio (fig. 5). 

There are now two separate sheets: the actual cover 
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sheet and the "substitute stiff ener sheet." Replace the 
two "by a single sheet located at the centroid of the two 
and with a thickness equal to the sum of the two thick- 
nesses. 

The fundamental formula (l) can now "be used to calcu- 
late the stiffness. 

The results obtained by this approximate method will 
approach those obtained with the more exact method as the 
amount of material in the stiffeners becomes smaller in 
relation to the amount of material in the cover sheet and 
also as the ratio !Li decreases. I*or the simple beams 



discussed in part II, the approximate results for tho 
stiffness were from 2 to 8 percent higher than thoso from 
the exact method. "" 

Nu merical examp le.- As a numerical example of the ge: 
eral case, the. stiffness calculation for' the schematic 
profile of figure 6 will be given. With 

A x = 89.3 A 3 = 250 A 3 = 125 

a l0 = 514.4 a 3D = 1875 a s0 = 5190 



= 166.7 



a 



S3 



= 333.3 



the equations become 



T = 2 (89.3 + 250 J s + 125 F 3 ) 




+ 166.7 (J a 






5190 P 3 + 333.3 (Jig - J a ) + 2U9 = 0 



The results are' 



I x = 0.001667 T 



? 3 . = 0.001254 T 
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F 3 = 0*000310 T 

Z&B - 0.01036 , 5? 

Substitution of these values in equations (6) and (8) 
gives the shear stresses indicated in figure 7, whero the 
direction of the stresses is also indicated, and the tor- 
sion constant J = 193*2 in. 4 . The torsion constant of 
the outer shell alone is J 0 = 113.6 in.*. 



II. COMPARISONS BETWEEN CALCULATION AilD EXPEKIUENT 



Te sts on box beam s.- Figure 8(a) gives the dimensions 
of three duralumin "box beams tested by the U.S. Army Air 
Corps (reference 6); figure 8(b) gives the dimensions of 
five stainless steel beams tested by the same agency (ref- 
erence 7). The test points for any given test run always 
fell very close to a straight line; individual test curves 
are therefore not reproduced in the present paper. 

The experimental values of the torsional stiffnesses 
are given in table I, which gives also the torsional stiff- 
nesses calculated by formulas (2) or by equations (5) and 
(8), The thickness t of— the vertical walls of tho stoel 
beams was replaced in these calculations by the effective 

thickness t e = 5. t because these walls formed diagonal- 

8 

tension fields at low torques. (See referonco 8.) 

The value of the shear modulus was taken as G = 4 X 10 

pounds per square inch for duralumin and Gr = 10 X 10 G for 
stainless steel. The shear modulus for stainless steol was 
taken lower than for ordinary steol in accordance with tho 
facts that the tension modulus of elasticity E for stain- 
loss steel is considerably Lower than that of ordinary 
steel and that the shear modulus of other highly alloyed 
steels has been found to be as low as 8 X 10 s ". 

Table I gives next the ratio of experimental stiffness 
to calculated stiffness. It will be Been that the average • 
ratio of all the beams is 0„90. 

In the group of stainless steel beams, beam 4 shows 
the smallest ratio, which may be attributed to the compli- 
cated section of the built-up stiffenors requiring many 
joints. 
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Duralumin "beam 1 has a ratio considerably in excess 
of unity. The test reports (references 6 and. 7) note that 
there must have "been some test irregularity, "because the 
experimental stiffness of beam 1 was higher than that of 
"beams 2 and 3, which "does not seem reasonable. 

Possible variations of elastic constants and of sheet 
thicknesses from the nominal values, however, may together 
account for perhaps 15 percent of the excess stiffness, 
which would reduce the experimental stiffness constant 
practically to the theoretical value. Since no direct er- 
ror could "be found that would make it necessary to elimi- 
nate the test rosult on beam 1 , it was included in the av- 
erage • 

A number of tests described in reference 6 and dis- 
cussed in reference 8 dealt with cases where the stiffness 
decreases very materially with increase of load. Tor the 
present purpose only the initial stiffnesses of thoso beams 
eoiao in question, and they are hot sufficiently well de- 
fined experimentally to give detailed numerical values; 
inspection of the results indicates that the ratio of ex- 
perimental to theoretical stiffness lies in general between 
0.9 and 1. 

It is suggested, .therefore, that T) = 0.90 be used 
as a general factor of effectiveness for box beams in tor- 
sion. ■ ' ■;; - - 

Ef fectiveness of stiffeners .- The outer cover of a 
shell is more effective in resisting twisting stresses 
than any material on the inside of the shell. It seemed 
of some interest, therefore, to calculate for the test 
beams the torsional stiffness of 'the outer cover alone, 
disregarding the stiffeners entirely, and to compare this 
stiffness ' &J Q with the calculated stiffness GJ, which 

includes the influence of tho stiffeners. The stiffness 
GrJ 0 . as well as tho ratio of G-J to G-j Q i3 givon in ta- 
ble I. Next in the table is* listed the reinforcement, 
expressed by the ratio of the area of the stiffeners to 
tho area of tho st if f enod ■ sheot . Piguro 9 shows the curve 
of increase in stiffness against reinforcement; it should 
be remembered that this curve is not generally applicable. 

A more logical way to represent the effect of the 
stiffeners, on the torsional stiffness is as follows: "im- 
agine the stiffeners removed; then increase the thickness 
of the sheet from which the stiffeners have been removod 
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■until — the torsional stiffness is the same as it was with 
tire st if f eners" present . Calculate the amount of. material 
that had to be added to the cover sheets; the ratio of 
this material to the actual stiffenor material gives a 
measure of the effectiveness of the stiffeners and is list- 
ed in the last row of tahle I. The efficiency of the hat- 
shape stiffener calculated on this hasis is only about 7 
percent. 

Te sts on. comp let e wings .- The experimental . wing de- 
scribed in reference 9 was chosen as the first example of 
a complete duralumin wing. The important characteristics 
of this wing are shown in figure 10. It consists of a 
central box with corrugated cover; but a smooth sheet 
forms a continuous cover over the nose, the central box, 
and the trailing-edge portion. 

The wing was first t-ested in the elastic range to de- 
termine the stiffness. Then it was subjected to the usual 
static tests, which were carried to. destruction in the 
high~ahgle-of~attack condition. The brxak was repaired, 
and the elastic stiffnesses were detimnined again in the 
following three conditions of the wing: (l) complete; (2) 
after removal of leading edge, leaving the central box and 
trailing edge; (3) after removal of- leading and trailing 
edges, leaving only the central box. 

Figure 11 shows the calculated and experimental .re- 
sults for the box alone. The agreement is very good ex- 
cept at the last — s"tart1ro~n .near the tip. Figure 12 shows 
the results for the combination of box and trailing edge; 
the calculated curve for the box alone is shown in dotted 
lines. It is apparent that the effectiveness of the trail- 
ing-edge pertion is very low, a. fact that can be easily ex- 
plained. The trailing-edge portion, buckles into a diag** 
onal-tension field. In order to realize the - theo ret ical 
stiffness of such a field, rigid flanges must be provided 
to take up the transverse component of the diagonal ten- 
sion. The trailing-edge strip, however, is so flexible 
that it furnishes practically no resistance tor these 
stresses and the cover sheet cannot develop any stiffness. 
Figure 12 shows that no large error is committed if the 
trailing-edge portion aft of the rear spar is entirely nsg- 
•lected;. this procedure . recommends itself also because it 
reduces very materially the labor of computing the stiff- 
ness of the complete wing, a two-cell box requiring very, 
much less computation than a three-doll box. The trailing 
edge was therefore neglected in calculating the stdffnoss 
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of the complete wing and figure 13 shows excellent agree- 
ment "between the calculated and the experimental twist, 

A point that deserves some mention is the calculation 
of the effective thickness of the webs. These webs had 
triangular lightening holes, making thejwebs, in effect, _ 
truss.es with double diagonals. The effective thickness of 
these webs was calculated "by formula (ll) of reference 10. 
This formula is valid for pin-jointed trusses; the trusses 
in question here have rigid joints, and the effect of 
those joints was taken into account by multiplying the 
calculated effective thickness by a factor 1.20 based on 
calculations and tests on similar trussed webs given in 
reference 3. 

Farther outboard, the webs had circular lightening 
holes. The effect of these holes on the shear stiffness 
of the webs was estimated from tests (reference 11). 

As the second example of a complete duralumin wing, 
the wing shown in figure 14 was chosen. It has a" central 
two-cell box covered with smooth sheet stiffened by closed- 
section stiffeners. The leading and trailing edges are 
fastened to the central box by means of piano hinges lead- 
ing them virtually open sections ineffective in torsion, 
at least at low loads. The torsionally effective part of 
the wing is indicated by the full lines in the plan and by 
tho cro sshatching in the sections. "" 

The interesting feature of this wing with respect to 
calculating the torsional stiffness is the large well for 
the retracting wheel. Obviously, only the box between the 
center and rear webs is ef f ective'j betwe en ribs 1 and 3. 
At rib 5, the full section from front web to rear web is 
effective. Between ribs 3 and 5, the front part of the 
box is cut up to form the recesses for the landing-gear 
struts. ITo detail drawings on which an estimate of the 
officiency could bo based were available for this portion; 
it was assumed, therefore, that the torsional stiffnoss 
varied linearly between ribs 3 and 5. 

The cut-out being close to the root, tho ioss in. tor- 
sional stiffness of the shell is partly balanced by "the 
bonding stiffness of the spars. The twist at rib 4 was 
therefore calculated by formula (10a) of reference 10, us- 
ing weighted averagos of moments of inertia and torsional 
stiffnoss between the root fittings and rib 4. It is oh- 
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viously not possible to calculate the curve of twist bo- 
tireen the root and rib 4, but this lack is of no practical 
imp'ortance , 

The test load was applied in four increments; the 
test points plotted in figure 15 are the averages of the 
four sets of readings. The outboard part of the shell 
buckled about in the middle of the t-est range; calcula- 
tions were therefore made for the nonbuckled state as well 
as for the buckled state. The agreement butween experi- 
ment and calculation is good in the outboard part; inboard 
there is' a rather large percentage error, but - this error 
is not large' in absolute magnitude; it may be partly due ■ 
to jig deflection. 

Stiffness under larg e torques *- The theoretical for- 
mulas discussed give the initial torsional stiffness under 
very low torques. IPo.r torques in the practical working 
range it will very often be necessary to apply corrections. 
One type of correction has been mentioned and used in the 
text: If the smooth covor buckles to form a diagonal ten- 
sion field, the actual sheet thickness t must be replaced 

by the effective thickness t e = § * ' T bis effective 

thickness maj-- be decreased much further at high loads (ref- 
erence 8), Corrugated or well-stiffened sheet can prob- 
ably "be assumed to lose not more than 5 or ip percent of 
its stiffness until failure occurs although there is no 
direct experimental evidence available to substantiate 
this claim. 

It— is apparent, then, that it is necessary before 
starting a stiffness calculation to. define the load range- 
in which it is to apply and to malce ; corresponding correc- 
tions if necessary to the formulas for initial stiffness. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Commit"tee for Aeronautics, 
Langley Field, Va., July 16, 1935. 
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TABLE I 



Experimental and Calculated Torsional Stiffnessee 





Duralumin beams 


Stainless steel beams 


Uumber of beam* 


i 


2 
fj 


3 
o 




u 






0 


Experimental stiffness 
GJ (10 6 lb. -in. 2 ) 


45.1 


41.5 


43.2 


20.2 


23.9 


29.8 


30.6 


20.2 


Calculated stiiineog 
GJ 


38 5 


53 a 


55.5* 


25 6 


25.4 




32 7 


23 5 


Ratio of fliperimentai. to 
calculated stiffness 


1.171 


.772 


.778 


.790 


.940 


.987 


.935 


,860 


Average rati^ (ti) 




.9C7 








.901 






Stiffnea? GJ 0 disregard- 
ing stiffeners 




23.2 




18.49 


17.4 


28.1 


32.3 


20.9 


Ratio of GJ 0 to calcu- 
lated GJ 




2.32 




1.38 


1.46 


1.08 


1.01 


1.13 


Ratio of stiffenar cross- 
sectional area to area of 
cover sheet 




4.10 




2.51 


2.52 


.52 


.52 


2.24 


Effectiveness factor for 
stiffeners (see text) 




0.41 




.27 


.37 


.37 


.06 


.09 


Duralumins G = 4 x 10 6 l"b. per eq..in. liax. 
Stainless Bteeli G = 10 x 1*> 6 lb. par aq.ln. 


toriiue: 


T = 5000 in.-lb. for beams 1 to 3 
5! = 4000 in.-lb. for beams 4 to 8. 



*By approximate method. 
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Figure 3.- Diagram for notation. 
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Figs. 4,5,6,7 




(a) (d) 
figure 4.- Box beam with stiffeners. 
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Figure 5.- Stiffener section. 




Figure 6.- Schematic profile. 
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Figure 7.- Stress sheet for schematic profile. 



Root 
fitting* 



Typloal seotion 8ootlon at rib 3 *wing 

All beam* symmetrical about center line 




5.12 



Beam 1. 
Web channels 0.040 
Longlt. oorr. sheet 0.040 
P-iTl/a D-3/4 



6.13 J 



Figure 14.- 

Wlng with 
wheal out-out. 
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Bean 4. 

Veb channels 0.006 
Btiffener 0.005 
Outer oover 0.005 
D- S/4 
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Beam 3. 
Veb obannels 0.040 
Longlt. oorr. sheet 0.040 
Outer coyer 0.012 (smooth) 
P« 3-1/3 D»3/4 



Beam 3. 
Veb obannels 0.040 
Longlt. oorr. sheet 0.040 
Outer oorer 0.013 
Oorrugated (P« 1-1/4, D-3/8) 
P- 3-1/3 D-3/4 




Beams 5*6 
Veb obannels 0.005 

Beam 5 6 

Stiffener 0.010 0.006 

Covering 0.005 0.013 
D- 3/8 



Figure B(a).- Duralumin box beams. (1-3) 



Beams 74 8 

Veb obannels 0.006 

Beam 7 8 

Stiffener 0.005 0.010 

OoTsrlng 0.013 0.006 
D- 3/8 



figure 8(b).- Stainless steel box beams. (4-8) 
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Figure 9.- Influence of stiffeners on torsional stiffness. 
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Figure 10.- Materiel Division 55-foot wing. 
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Figure 11.- Twist of 'box alone. 
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Figure 13.- Twist of combination box and trailing edge. 
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Figure 13.- Twist of complete wing. 
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